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In the present work, we investigate the intrinsic relation between quantum fidelity susceptibility (QFS) and the 
dynamical structure factor. We give a concise proof of the QFS beyond the perturbation theory. With the QFS in 
the Lehmann representation, we point out that the QFS actually the negative-two-power moment of dynamical 
structure factor, and illuminate the inherent relation between physical quantities in the linear response theory. 
Moreover, we discuss the generalized fidelity susceptibility (GFS) of any quantum relevant operator, that may 
not be coupled to the driving parameter, present similar scaling behaviors. Finally, we demonstrate that the QFS 
can not capture the fourth-order quantum phase transition in a spin-1/2 anisotropic XY chain in the transverse 
alternating field, while a lower-order GFS can seize the criticalities. 
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I. INTRODUCTION 


Recently the concept of quantum fidelity susceptibility 
(QFS) has been recognized as a versatile indicator in identi¬ 
fying quantum critical points (QCPs) III. Plentiful literatures 
reveal that the QFS exhibits singularity at a great variety of 
critical points Nevertheless, there still remains vague 

issues in understanding the QFS. For a many-body Hamilto¬ 
nian H (A) = Hq + XHj, when the driving parameter A varies 
across a transition point Ac, the noncommutability between 
Hq and Hj leads to a nonanalytic behavior in the ground state 
'Po and a reconstruction of low-energy spectra Ei. Originated 
from quantum information theory, quantum fidelity is defined 
as the overlap of wavefunctions at two points in parameter 
space: 

.F(Ao,Ai) = A/(^o(Xi)I^o(Aij)(^o(Aij|^o(Ao)y. (LI) 


It was discovered that quantum criticalities promote the de¬ 
cay of fidelity d. Such sensitivity to the quantum criticality 
was further amplified by the associated changing rate. In this 
respect, the QFS is defined as the leading-order of the Taylor 
expansion of the overlap function F{X,X + SX) d, 


XF 


-2\nF{X,X + SX) 

A™ (JA)2 


( 1 . 2 ) 


As A crosses QCPs adiabatically, xf shows a peak singularity 
and exhibit universal scaling law. Such universal properties of 
QCPs in equilibrium prompted a rapidly growing theoretical 
interest to the quench dynamics. 

When the tuning parameter A is time dependent, say, A(t) ^ 
vf^, where v is the adiabatic parameter that controls the prox¬ 
imity to the instantaneous ground state, the energy gap A at 
the critical point Ac tends to zero when the size of the system 
N ^ oo and it is almost impossible to pass the critical point 
at a finite speed v without exciting the system. The nonana¬ 
lytic scalings can be revealed by singularities in the general¬ 
ized adiabatic susceptibility X 2 r-i- 2 (A) of order 2r + 2|il: 


^ |(^KA)|g/|^o(A))p 

^ [Ei{X)-EQ{XWr+2^ 


(1.3) 


where ip; and Ei (I ^ 0) correspond to the excited states and 
energies respectively. The general theory for crossing the crit¬ 
ical point at a finite rate is given by Kibble-Zurek mechanism 
mi 113]. In the opposite limit, i.e., u—/ 0, Ref. fllll firstly 
discussed the adiabatic evolution of ground state in the frame 
of perturbation theory, and obtained 


^ |(vPz(A)|gz|vPo(A))p 


(1.4) 


This exactly corresponds to the generalized adiabatic suscep¬ 
tibility [Ea.( IL3b 1 of order 2. Note that the second derivative 
of the ground-state energy is of order 1. 

Due to the arbitrariness of relevance of the driving Hamilto¬ 
nian Hi under the renormalization group transformation, xf 
increases as the system size grows and manifests distinct scal¬ 
ing behavior. The summation in Eq. (IL41 i contributes to an 
extensive scaling of xf in the off-critical region. Therefore, 
QFS per site xf/N appears to be a well-defined value, where 
N = L‘^ is the number of sites and d stands for system di¬ 
mensionality. Instead, the QFS exhibits much stronger depen¬ 
dence on system size across the critical point than in noncrit- 
ical region, showing that a singularity emerges in the summa¬ 
tion of Eq. CT. This implies an abrupt change in the ground 
state of the system at the QCP in the thermodynamic limit. 
For example, when Hi is charge-or spin current operator J or 
the polarization operator P, it is shown that in gapless one¬ 
dimensional (ID) systems with open boundary conditions the 
leading terms in the L dependence are given hy xj oc KL^ 
and Xp oc KL^/u^, where K is the Luttinger liquid param¬ 
eter, u is the characteristic sound velocity, and the numerical 
prefactors are universal S. Following standard arguments 
in scaling analysis of a second-order quantum phase transi¬ 
tion (QPT), one obtains that the QFS per site in d dimensions 
in most cases {vd > 2) scales as lUAllall^ 

Xp/A^~L'/"-‘'/(|A-A,|L'/"), (1.5) 


where v is the critical exponent of the correlation length and 
/( .) is a regular scaling function. As for the first-order QPT, 
the QFS per site scales exponentially Cl, i.e., 

Xf/N ^ g{N)e^^, 


( 1 . 6 ) 
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where ^ is a size-independent constant, and g{N) is a poly¬ 
nomial function of N. 

In addition to an increasing interest in the finite-size scaling 
analysis, a few researches addressing that the QFS cannot de¬ 
tect the QPTs of infinite order, have joined in the subject of ac¬ 
tive study fllUlIi. This conclusion conflicts with the find¬ 
ings on the ID Luttinger model and the XXZ chain 

by virtue of the density-matrix renormalization-group tech- 
nique ll^ and the real-space quantum renormalization group 
but supports initial opinion proposed by You et al. El. 
However, recently G. Sun et al. argued that using the QFS as 
a sensor to detect Berezinskii-Kosterlitz-Thouless-type tran¬ 
sitions is uncontrollable 1^ . Therefore, the issue on fidelity 
susceptibility in describing high-order continuous phase tran¬ 
sitions seems still controversial. 

Contrary to the QFS, the thermal-state fidelity is deter¬ 
mined as a descendant of Uhlmann fidelity ll2l 

Ft{pq, Pi) = tr^ y/^Poy/^, (1-7) 


d\H’s within the adiabatic perturbation theory. To this end, 
the susceptibilities can be also expressed through the imagi¬ 
nary time correlation function of the perturbation Hj (r) using 
the relation El l20l] : 

nOO 

XF= / r[(T-o|i7/(T)i7/(0)|T'o)-(«'o|i7/|«'o)']dT, 

Jo 

( 2 . 1 ) 

where iT/(T)=exp(iFT)iJ/ exp(—iTr). In this spirit, the QF S 
was extended straightforwardly to finite temperature Ola 1^ ; 

Xf(/ 3)= / r[(T-o|iT/(T)iJ,|T'o)-(^'o|i7/|^o)']dT, 

Jo 

( 2 . 2 ) 

in which /3 is the inverse temperature. 

In the Lehmann representation lEg. lll.Tl il. the hopping ma¬ 
trix between ground state and arbitrary excitations is written 
as 


where pi (i=0,1) are density operators. A simple formula is 
achieved by considering temperature driven phase transitions 


Ft{Po,IJi) 


Z[(/3o+/3i)/2] 
v/X(/3o)X(/3i) ■ 


( 1 . 8 ) 


Then there is a direct connection between the thermal fi¬ 
delity susceptibility and a purely thermodynamic quantity 
El , namely, specific heat. 


(T'/|iT/|T'o) = = {Eo - Ei){'i'i\dx-^o)-(2.3) 

Putting the above relation into right hand of Ea. (ll.41 i and mak¬ 
ing use of completeness relation = 1, we obtain 

E = (dx^oldx^o) - |(5A4/o|'ko)p.(2.4) 

Remarkably, the QFS was now devised to seize the criticality 
in the perspective of the Riemannian metric tensor formEl 


XT = 


4/3- 


2 ■ 


(1.9) 


Eq. (11.9b can not fall to the QES in zero-temperature limit, 
P ^ oo. Meanwhile, the QES for degenerate ground state is 
ill defined. 

In the present work, we analyse the immanent relation be¬ 
tween the QFS and the dynamical structure factor. We give 
a concise proof of Eg. ( 11.4b beyond the perturbation theory 
in SecHI] With the QES in such Lehmann representation, we 
point out that the QES actually the negative-two-power mo¬ 
ment of dynamical structure factor in SecHIIl and illuminate 
the inherent relation between physical quantities in the linear 
response theory. In Sec lIVI we discuss the critical behavior of 
the generalized fidelity susceptibility (GES) in ID transverse 
Ising model. We also show that the QFS fails in detecting 
the QPT in ID anisotropic XY model under a transverse al¬ 
ternating field, while a lower order GFS Ad(_ 3 ) can probe the 
fourth-order QPT therein. The paper is summarized in Sec. 

lY] 


II. FIDELITY SUSCEPTIBILITY IN NONPERTURBATIVE 
FORM 


The universal infrared divergence comes from the low- 
frequency contribution in Eq. ( IL4b . This prompts us to con¬ 
sider the dynamical response of the system to the interaction 


Xf = (Oa^'oIOa'I'o) — |(i9a4'o|^o)|^. (2.5) 

The implications of the equivalence between Eo. (l2.5b and 
Eo. (IL4b are three fold. Eirst, it shows that the QFS is well de¬ 
fined as a consequence of orthogonality of energy levels, with¬ 
out prior knowledge of fidelity. Second, the usage of Eo. (l2.5b 
is nonperturbative, so the adiabatical perturbation theory is not 
indispensable. Last but not the least, the description of Rie¬ 
mannian manifold abandons the requirement that the parame¬ 
ter A is coupled with the driving Hamiltonian. In other words, 
we can investigate the fidelity susceptibility of any quantum 
operator O, that may not appear in the Hamiltonian explic¬ 
itly. Such tentative investigation has been demonstrated in 

Ref.ii. 


III. DYNAMIC STRUCTURE FACTOR 


A most common experimental characterization in con¬ 
densed matter physics is inelastic neutron scattering, which 
directly measures the energy and momentum dependence 
of spin-spin correlations as described by the momentum- 
resolved dynamical structure factor. For a quantum operator 
O, the momentum integrated dynamic structure factor S'jw) is 
defined as ll^ 

S{d,L0) = Y, |('I^(|0|vI/o)P<5(u; - wjo). (3.1) 
//o 
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Here w/o = Ei — Eq denotes the corresponding excitation 
energies of Zth excited state. S{uj) contains rich information 
of a many-body system. For instance, the dynamic response 
function in linear response theory can be expressed asyH] 


pOO 

X(0,uj) = / duj' 

Jo 


S(d,u;') 
w + w' + i0+ 


S(6^,u;') 

LU — LU' + i0+ 


(3.2) 


which can be also linked with experimental signal by apply¬ 
ing external perturbing field through Kubo formula. If O is 
Hermitian or the system is time-reversal invariant, as is sup¬ 
posed in the following, |(^'/|0|'I'o)P = |('I'z|0^|'I'o)P> then 
Eg. (13.21) takes on a simplified form: 


where the summation runs over both I and n excluding they 
are equal. Here partition function Z = g-/3s„ 

(jjin = El — En denotes the energy difference between ex¬ 
cited states 1^) and \n). In the zero-temperature limit, \n) is 
restricted to ground state |0), and the dynamic form factor co¬ 
incides for positive ui, since only processes by which energy 
is transferred to the system are allowed. At finite temperature, 
the corresponding GFS of the dynamic form factor 


■^(fe)((0)) = / dcouj’^Srid^u;) 


n^l 


(3.10) 


pCO 

x{O, 0 j) = / 2(jj'dw' 
Jo 


S{d,u;') 


uj'^ - (w -I- i0+)2 


(3.3) 


The linear response function can be related to the moments 
Al(fc) of the dynamic form factor, defined by 


pOO _ 

M(k)(d)= / doju;^S(d,u;) = yojto 
Jo , 


(3.4) 


where k is an integer. The moments of different orders give 
rise to various sum rules and have explicit physical meanings. 
For example, zero-order moment, also called as static struc¬ 
ture factor, quantifies the fluctuation at zero temperature, 

M(o)(d) = (^'olO^l^o) - |(^o|0|^o)|", (3.5) 

which is the leading response in linear quench dynamics 1^ . 
Similarly, negative-one-order moment corresponds to static 
response function 

M(_i)(d) = ^x(O,uj = 0). (3.6) 

From Eq.( l2.5l ). negative-two-order moment is the QFS with 
the relation 


Note that these GFSs, contrary to the T = 0 case, get contri¬ 
bution from negative uj as well, since transferring energy both 
from the particle to the system, or from the system to the rest 
particle are both admissible at T > 0. 

Much attention has been paid to the critical behavior of cor¬ 
relation function in the QPT In the low-frequency limit, the 
dominated contribution comes from static response function 
X. Provided x shows a clear singularity, such anomaly can 
act as a signature of criticality. While for some phase tran¬ 
sitions, X does not manifest singular trend, and then Ad(_ 2 ), 
equivalently the QFS, can be used as a supplementary means 
of identifying the phase transitions. The QFS was proven to 
successfully characterize the critical behavior of QPTs in a 
many-body system H. 

Common speaking, the reconstruction of low-energy exci¬ 
tations mainly contributes to the negative-order GFS while the 
change in high-energy part is crucial to positive-order GFS. 
When quantum operator O acting on ground state has nonzero 
components with low-energy states, the singularity is more 
pronounced for larger magnitude of |fc|. In fact, for arbitrary 
k, we have 


^min 


< 


M(k+i) 


(3.11) 


3W(-2)(0) = xf(0). (3.7) 

Eg. (13.71) bridges the QFS and dynamic structure factor, and 
enriches the techniques to obtain the QFS. In other words, Eq. 
( 13.4b can be seen as a generalization of the QFS, and Ad(fe) 
is named as fc-order GFS. Although S{0,uj) is not easy to 
retrieve, the dynamic structure factor x(w) can be solved in 
terms of a series of Feynman diagrams, and thus the dynami¬ 
cal structure factor can be gotten from Eg. (13.2b : 

S(0,„) = bw(^, („>0), (3.8) 

TT 

As far as S{0, w) is known, finite-order GFSs can be gained. 

For a quantum operator O at equilibrium at a finite temper¬ 
ature T, the momentum integrated dynamic structure factor 
St{oj) is defined as 


where Wmin is the lowest excitation. The equality only holds 
for the existence of fl'min and (4'min|0|T'o) = 1. It is clear 
that the lower power of negative GFS displays more diver¬ 
gence provided gap closes at criticality. 

IV. APPLICATION AND RESULTS 

To see such generalization in Eq. ( 13.4b plays an essential 
roles in various types of QPTs, we calculate the GFSs in a 
few examples in the following. There has been intense interest 
in studying ID spin systems such as transverse Ising model, 
which serve as a general prototype for quantum magnetism 
and for the understanding of QPTs. All the concepts about 
both equilibrium and non-equilibrium QPTs have been tested 
on this model. The ID spin-1/2 transverse Ising model is well 
known for its solvability and paradigm, given by 


^ ^__ p PJJ'n 

St{0,u;) = ^ (3.9) 

l^n 


N 


i=i 


j ^3 + ^ 




(4.1) 
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where h is the strength of transverse field and N is the num¬ 
ber of spins (assumed here to be even). The Hamiltonian ( 14.1b 
respects the Kramers-Wannier duality symmetry, reflected in 
that the site in the original chain has one-to-one correspon- 


dence with the bond of the 

dual chain and vice 

versa by a 

non-local mapping: 



7 = 


(4.2) 

7 = 

n-t- 

(4.3) 


k<j 



The operators defined in Eqs. dOi-dOi satisfy the same set 
of commutation relation as the Pauli operators, i.e., they com¬ 
mute on different sites and anti-commute on the same site. 
The Jordan-Wigner transformation maps explicitly between 
pseudospin operators and spinless fermion operators by 


cr = exp 


cr- = exp 


j-i 

•E' 

i=l 

f -1 


3-1 


-m^c\ci 


2 = 1 


j -1 


= n 

2 = 1 




l-2c]c,. 


(4.4) 

(4.5) 

(4.6) 


After the fermionization, Fourier transformation can be pro¬ 
ceeded consequently. 


1 

7 ^ 



q 


(4.7) 


where Wg = Such that the Hamiltonian 

becomes a quasi-free fermion system 

= (4.12) 

q 

Obviously the ground state energy is obtained, 

Eo = -Y,^g. (4.13) 

There is an energy gap in the excitation spectrum of the sys¬ 
tem which goes to zero at g = 0 for h = 1, i.e., 

A{h) = El - Eo = 2\1 - h\, (4.14) 


indicating the divergence of the correlation length and a QPT 
at /ic = 1 from an ordered state ((cr®) = 1 ) to a disordered 
state ((cr®) = 0 ), which is the consequence of ferromagnetic- 
paramagnetic duality. To investigate the held induced criti¬ 
cality, a natural choice is z-component magnetization 0=a^. 
The dynamical structure factor of spin magnetization yields 

S{a\uj) = Y,^-^S{u3-2u3,). (4.15) 

g ^<1 


The static response function takes on 

Jo ^ 


= E 


sin^ q 


. (4.16) 


Away from critical point ft, 7 ^ 1, we have a linear scaling as 


with the discrete momenta q, which correspond to the reduced 
Brillouin zone and are given by 


nn 


9 =—, n = -(iV-l),-(iV-3),...,(iV-l). (4.8) 

Hereupon the Bogoliubov-de Gennes (BdG) Hamiltonian is 
obtained as 


H = - 


1^(2 cos q - 2h)c\cq + i sin + c-qCq) (4.9) 


where K[a\ and E[q\ correspond to the first and second kind 
complete elliptic integrals of a = 4ft/(l + h)^. Around the 
critical point ft = 1 , we have 




N 

x(o■^0)U^l 


In I ft — 11, 
A^lnAT. 


(4.18) 


Bogoliubov transformation can be adopted to diagonalize the 
quadratic Hamiltonian via 

Cq = COs9qbq + ism9qb^_^, 
cj = cos9qb^q — ism9qb-q, 


C-q = cosftgft_g — isinftgftj, 
(J_q = COSftgftLn + isinftgftg, 


(4.10) 


where bq and are fermionic operator satisfying the same 
anticommutation relation as Cg and cj. The coefficients are 


determined by 


cos 29q = 


sin 29 q = 


cos q — h 

LUq 

sing 

U3„ 


(4.11) 


The QFS is simultaneously given by 

S{a^,u}) 




dco- 


= E 


sin^ q 

g 9 


(4.19) 


The result ( 14.19b is identical with that from Fa . (12.5b directly. 
A closed form of the QFS was revealed in Ref. 11^ : 


Xf( 6 -^) = 




'-N 


N 


- ft^ 


16ft2(ft^-bl)2 16ft2 (ft^ -b l)(ft2 - !)■ 

(4.20) 


Away from critical point ft 7 ^ 1, we have 

Xf ( ct ^,0 < | ft | < 1 ) 1 


N 

XF(<7Mft| > 1 ) 
N 


16(1 -ft2)’ 
1 


16 ft 2 (ft 2 — 1 ) ’ 


(4.21) 

(4.22) 

























5 


and around critical point. 


_,2 


1 ) = ^ 


cot^(g/2) 1 

16 ^q’ 


(4.23) 


where the lower bound gives the leading contribution in the 
cotangent series. Thus v = 1 can be retrieved by Ea. (ll.51 l. 
Particularly, Kramers-Wannier duality of the Ising model 
grantees the relation: 

M = (^)^XF(<5•^ i). (4.24) 

Figure [T] shows that apparently the QFS exhibits more pro- 



h 


FIG. 1: (Color online) The magnetic susceptibility 14.16b and fidelity 
susceptibility 14.19b as a function of magnetic field h on N = 100 
transverse Ising chain. 

nounced singularity than magnetic susceptibility. Remark¬ 
ably, Ref. 11^ disclosed a potential measurement of the QFS 
through the zero-momentum dynamical structure factor. 

There are other choices of O, e.g., momentum space Majo- 
rana operator Omf = The dynamical structure 

factor of Majorana operator yields 

^((^mf.w) = ^(5(a; - w,). (4.25) 

q 


Accordingly, 


(jJt UJt 


The asymptotic behavior appears that 
Xf(Omf) 1 


N \l-h?y 

Xf(Omf, —>■ 1) ~ iV^. 




(4.27) 


ward to understand it is impossible to find an observable op¬ 
erator O, which is not related to driving parameter A and ex¬ 
hibits more singular QFS than that of (9 mf, because the de¬ 
nominator of Fa. (ll.4b is at most inverse square of gap ( 14.14b . 
that is. 


Xf( 0) < ^ E l('f'i(A)|0|vko(A))p. (4.28) 





h 


FIG. 2: (Color online) The GFSs Af (_ 3 ) of and Majorana oper¬ 
ator (5mf as a function of magnetic field h on N = 100 transverse 
Ising chain. Al(_ 3 ) (Omf) has been divided by 16 to make it com¬ 
parable with Af(_ 3 )((T^). 

Provided xf is not sensitive enough to observe a QPT, 
GFSs Ai(k) with more negative k (k < —2) is needed, such 
that 


q q 

7\d(_3)((9MF) = E (4.30) 

q 9 

Away from critical point, we have Ad(_3)((5) ~ N, in con¬ 
trast to Al(_ 3 )((i)) 3%. 7V^ sitting on the critical point. The 
GFSs A'l(fc) with more negative k takes on more singular 
peaks, shown in FiglJ] To validate the generalization to lower 
order of GFS is significant, we resort to another less investi¬ 
gated model. 

The ID spin-1/2 anisotropic XY model in a transverse al¬ 
ternating field is defined by 


with 


hj = h- {-yS, (4.32) 


A natural question is that whether we can find operator O, 
whose scaling exhibits stronger singularity. It is straightfor- 


where 7, h and S are anisotropy in the XY plane, the uniform 
and alternating magnetic field strengths, respectively. Here 
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periodic boundary conditions are assumed. Using Jordan- 
Wigner transformation in Eqs. ( l4.4H4.6l l, the model can be 
diagonalized as: 

N 

H = - X]{7(c+ C ++1 - CjCj+i) + (c+ Cj+i - CjC+^i) 

3 

- [h-{-y5]{l-2c+c,)}. (4.33) 

We introduce the discrete Fourier transformation of two kinds 
of fermionic operators in this period-two system: 

where N' = N/2 and the discrete momenta are given as fol¬ 
lows: 

71 TT 

k = —, n = -{N'- 3),..., {N'- 1).(4.35) 


In terms of vector operators Pj' = (a|,, a_fc, 6^, 6_fc), we 
rewrite it in the BdG form, 


H = 

k 


(4.36) 


where 


nk 


^ {h + S) 0 
0 -{h + 6) 

7(1-6" = '=) (i-r 6 -'=) 

\ 2 2 


2 

7 ( 1 - 6 ='=) 

{h-5) 

0 


7 ( 1 - 6 ='=) 

2 

(i+£^ 


2 


\ 

(U37) 


0 

-{h-5)l 


The diagonalization of the BdG Hamiltonian (14.371) has un¬ 
dergone an exhaustive study. More precisely, the diagonalized 
form is achieved by a four-dimensional Bogoliubov transfor¬ 
mation which connects {uk, O-fe, &fe, &-fe} with two kind of 
quasiparticles. 


7fc.i 

7lfc,i 


Ik,2 

\ 7lfc,2 / 


/ ak \ 

alk 

bk 

V blk 


(4.38) 


where the rows of Uk are eigenvectors of the BdG equations, 

[H,llj=ek,111,, ( 4 . 39 ) 

[H,ll2] = Sk,2ll,. (4.40) 

Here ek,i and ek ,2 are elementary excitations. To this end, 

^ = EE2£fc."(7d7fe.n-2), (4.41) 


where 

Sk,l = 2 




— 2 cos 2 ^-I-(52(h2-I -72 sin^ 


1/2 


(4.42) 


£k,2 = 2 


h^+S^ + cos^’l+^^si^^^ 


+ 2 \I h'^ cos2 -I- + 72 sin^ —) 


1/2 


(4.43) 


The ground state energy is thus obtained, 

=^0 = - E (^'=-4 + efc. 2 ) ■ (4.44) 


Since QPTs are caused by nonanalytical behavior of ground- 
state energy, QCPs correspond to zeros of £k,i. Quantum 
phase boundaries are determined by the equations: 

+ j 2 ^^ 2 _^^ 2 ^ (4 45^ 

Apart from a second-order broken-symmetry QPT at hc,i = 
iV(52-1-1 [see Fig.(l^)], it was discovered that a fourth-order 
QCP occurs at (h —>■ 0,(5 = ± 7 ) and [h = ±1,(5 —?> 0) ll^ 
[see Fig.(12i)]. A critical behaviour is characterized by a set of 
exponents which determine peculiarities of different ground- 
state quantities in the vicinity of the critical field. The corre¬ 
lation length ^ oc |h — and the transverse static sus¬ 

ceptibility diverges as (x |h — hc|"“, and the energy gap 
vanishes as A oc \h — , where a and z are the critical 

exponent of specific heat and the dynamic critical exponent. 
The transverse static susceptibility, i.e., the second-order dif¬ 
ferential of the ground-state energy, = d'^Eo/dh'^ diverges 
at hc,i = ±-\/(5^ ± 1. FigurelUreveals that oc |/i — hc,i\~°‘ 
with Of = 0 , that is to say, only a logarithmic divergence in 
the dependence of h jj^ . The elementary excitation energy 
occurs at k = 0 which permits us to get for the energy gap 
in the vicinity of the critical points, A oc |/i — imply¬ 
ing izz = 1. When 7 > 0, generic QCPs belong to the two- 
dimensional Ising universality class with critical exponents iz 
= 1, z =1. Nevertheless, we find that around he ,2 = 0, 

(x \h — hc, 2 \~°‘ with a = —2, i.e., x^ does not diverge at 
QCP at hc ,2 and d'^x^l^h'^ exhibits a logarithmic divergence 
at hc, 2 , implying a fourth-order QPT. The energy gap occur¬ 
ring at fc = TT in the vicinity of the critical points behaves as 
A cx |hp, implying izz = 2. The critical exponents v = 2, 
2 ; = 1 are identified via careful analysis, corresponding to an 
alternating universality class. The obtained critical exponents 
confirm one of the scaling relations in QPT: 2 — o; = {d + z)v. 


By Fq. (13.41) . the GFSs of different orders are showcased 
in Fig. m We find that the critical points hc,i have been al¬ 
ready captured by the singularpeakof while only 

Ad(_ 3 )(CT^) can detect the anomaly around he ,2 with a singu¬ 
lar peak. We can observe that a similar behavior takes place 
for the GFSs of Omf; see Fig. 0 We can speculate more 
pronounced peaks of Ad(fe) (O) will emerge for k < —3. 
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FIG. 3: (Color online) Derivatives of ground-state energy with re¬ 
spect to h in the ID spin-1/2 anisotropic XY model under a trans¬ 
verse alternating field, (a) -d^Eo/dh^\ (b) -d'^Eo/dh*. Parameters 
are as follows: 7 = 0.5, <5 = 0.5, N = 8000. 


FIG. 5: (Color online) The GFSs of Omf with different orders in 
the ID spin-1/2 anisotropic XY model under a transverse alternating 
field. Parameters are as follows: 7 = 0.5, 5 — 0.5. 



0 

h 


FIG. 4: (Color online) The GFSs of with different orders in the ID 
spin-1/2 anisotropic XY model under a transverse alternating field. 
Parameters are as follows: 7 = 0.5, 5 = 0.5. 


V. CONCLUSIONS AND DISCUSSION 

In this paper, we proved the equivalence of two forms of 
the QFS between Eg. (11.41) and Eg. (12.51) . We show in a sense 
that the QES is neither a pure descendant of the concept of fi¬ 
delity, nor a response of the perturbed Hamiltonian. Instead, it 
connects the critical properties with many relevant operators. 
We discuss the QES of z-component magnetization and 
momentum space Majorana operator (9 mf in the ID trans¬ 
verse Ising model, which presents similar critical behaviors. 
The QES can be generalized to moments Ai(^k) of the dy¬ 
namic form factor, which are experimentally accessible. In 
the vicinity of a critical point, a lower-order GES A4(fc_i) 
{k < —2) will exhibit stronger singular behavior than that 
of Considering that a singularity in the second deriva¬ 

tive of energy seizes a second-order QPT, it is supposed that 
of relevant operators is requisite to detect the crit¬ 
icality for a fc-order phase transition. We can observe that 
A4(_i), namely, magnetic susceptibility has detected the 


second-order QPT in ID transverse Ising model, and lower- 
order GESs display much stronger peaks. Nevertheless, for 
the fourth-order QPT occurring in the ID spin-1/2 anisotropic 
XY model under a transverse alternating field, we find that 
Ad(_ 2 ) does not show any anomaly in the vicinity of QCP, 
while A4(_3) or lower-order GESs can capture the singular¬ 
ity. This finding infers that the QES may not predicate the the 
existence of the Berezinskii-Kosterlitz-Thouless QCP, which 
affirms our argument in Ref. d from another perspective. 

In the formulas (13.4b and (13.10b . the structure change is en¬ 
coded in the transition matrix elements |(T';|(5|'l'o)P of an 
observable O between the ground and excited states. The 
scope of the GPS can be reaily applied to the related solv¬ 
able models like the transverse XY chain dimerized XX 
model IH], the BCS Hamiltonian 14^ . the Harper chain EH], 
the Kitaev model S and compass chain 15 - 1^ . Also, 
the concept of moment is easily extended to finite temper¬ 
ature Besides in those solvable models or very small 

systems, the matrix elements can also be potentially updated 
and evaluated via the prevailing numerical techniques, like 
the Monte Carlo methods and transfer-matrix renormalization 
group technique. More importantly, in experiment the fre¬ 
quency resolved structure factor is easily measured via the 
well known techniques, such as neutron scattering, angle re¬ 
solved photoemission spectroscopy (ARPES) techniques, etc. 
Por instance, the ID quantum transverse Ising model was re¬ 
alized in the insulating ferromagnet CoNb 206 under strong 
transverse magnetic fields and the spin dynamics was mea¬ 
sured by neutron scattering El]. By Eq. dM]), the GPS can 
be derived from the knowledge of the dynamic structure fac¬ 
tor. 
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